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Abstract 

It has been estimated that the entire Earth generates heat corresponding to about 40 TW (equiv- 
alent to 10,000 nuclear power plants) which is considered to originate mainly from the radioactive 

OO . 

, decay of elements like U, Th and K, deposited in the crust and mantle of the Earth. Radioactivity 

cn 

\ of these elements produce not only heat but also antineutrinos (called geo-antineutrinos) which can 

Q^' be observed by terrestrial detectors. We investigate the possibility of discriminating among Earth 

composition models predicting different total radiogenic heat generation, by observing such geo- 



antineutrinos at Kamioka and Gran Sasso, assuming KamLAND and Borexino (type) detectors, 
respectively, at these places. By simulating the future geo-antineutrino data as well as reactor 
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■ antineutrino background contributions, we try to establish to which extent we can discriminate 

among Earth composition models for given exposures (in units of kt- yr) at these two sites on our 
planet. We use also information on neutrino mixing parameters coming from solar neutrino data as 
well as KamLAND reactor antineutrino data, in order to estimate the number of geo-antineutrino 
induced events. 
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I. INTRODUCTION 



There is much about the Earth's heat engine that is unknown. The main datum is the 
surface heat flow, as most of the Earth is hidden from view. The total heat flux is presently 
estimated to be about 40 TW, which, however, suffers from uncertainties due to the size of 
its local variations and inaccessibility of much of Earth's surface. 

We know that there are radioactive isotopes in the Earth that can produce heat through 
decay. Although the rate of heat generation by decay of unstable radioactive nuclides is 
tiny, if it is integrated over the entire volume of the Earth, the total heat flux becomes 
huge. Radiogenic heat, evidently, must be an important source of internal heat production. 
However, questions as how much of the Earth's heat generation is from radiogenic origin, 
how much from residual heat remaining from the formation of the Earth or from the release 
of gravitational energy as the Earth contracts, have not yet been answered in a satisfactory 
way. 

In the context of modeling the thermal state and thermal history of the Earth it is 
important to know the spcciflc heat production of the chief heat-producing radionuclides 
such as '"'K, ^^^U, ^^^Th and ^''Rb, encountered on surface layers and supposed to exist also 
in interior layers. The heat that drives mechanical motion in the mantle presumably comes 
mostly from radioactivity. In the radioactive decay process, a portion of the mass of each 
decaying nuclide is converted to energy. Most of this is the kinetic energy of emitted a and 
P particles or of 7 rays, which is fully absorbed by rocks within the Earth and converted to 
heat. 

For (3 decays, however, part of the energy is carried away by emitted neutrinos and 
antineutrinos. Clearly, measuring the neutrino and antineutrino fluxes from the Earth can 
provide a unique way to access information on the internal structure and dynamics of our 
planet [1-9]. 

This can be done, at least for electron- antineutrinos {Vf.) coming from ^^^U and ^32^]^ 
decays, the former producing 6 and the latter 4 z/g per decay chain within the energy reach of 
current and near future neutrino detectors. Antineutrinos from these elements have higher 
energies than those that come from and ^''Rb, and are detectable at liquid scintillator 
detectors such as KamLAND [10] and/or Borexino [11] through the inverse /3-decay reaction. 
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p + Uf, n + where the detection threshold energy is Ej^^ = 1.7 MeV. 

Different Earth composition models predict different total amount of U and Th in the 
mantle, which lead to different total heat flux from radiogenic origin. Moreover, the fact that 
concentration of such elements is much larger (a factor of ~ 20 times) in the continental crust 
(with typical thickness ~ 35 km) than the oceanic one (with typical thickness ~ 7 km) and 
continents and ocean are not uniformly distributed over the Earth make the geo-antineutrino 
flux different from place to place. 

Quite recently, in Ref. [7], geo-antineutrino fluxes from various Earth composition mod- 
els have been estimated (before the first KamLAND results were reported) and then in 
Ref. [8] , the first KamLAND data, which contained ~ 9 possible geo-antineutrino candidate 
events [10], have been analyzed, and it was concluded that practically all Earth composi- 
tion models are consistent with the current data. Moreover, the accumulated data so far at 
KamLAND is not enough to establish the presence of geo-antineutrinos [10] and we must 
wait for future data. 

In this work, we will try to go beyond the above mentioned papers, by investigating to 
which extent hquid scintillator (type) detectors such as KamLAND and/or Borexino can 
be used to help in discriminating among different geophysical models of heat production 
by measuring antineutrinos in the energy range 1.7 < Ep^/MeV < 3.4 produced inside the 
Earth. We show how much one can improve the quantitative understanding of the radiogenic 
contribution to terrestrial heat in about a decade of operation. 

The organization of this paper is as follows. In Sec. II, we briefiy describe a few models for 
the Earth radioactive element composition, in Sec. Ill, we present our analysis method and 
in Sec. IV, we describe our results. Finally, Sec. V is devoted to discussion and conclusion. 

II. EARTH AS A ANTINEUTRINO SOURCE 

Antineutrinos are produced inside the Earth in radioactive (3 decays mainly of ^°K, ^^^U, 
^^^Th and Rb. These elements are classified as lithophile elements in geophysics and 
considered to be accumulated more in the Earth's crust. The abundance of these isotopes, 
although of prime geophysical importance, is only known at or near the surface of the Earth. 
Among these geo-antineutrinos, the one come from ^^^U and ^^'^Th have higher energies. The 
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maximal neutrino energy from tiie former and the latter are, respectively, E?^^^ = 3.26 MeV 
and E^^^ = 2.25 MeV, which are above the threshold of existing anti-neutrino detector such 
as KamLAND [10] and Borexino [11]. 

Here we consider, as our references, four different models which predict the distributions 
of ^^^U and ^^'^Th in the Earth. We do not consider neutrinos coming from ^'^K and ^^Rb 
because energies of these neutrinos are below the threshold of detectors we consider in 
this work. For all models, we have considered three regions, namely, continental crust, 
oceanic crust and mantle, with different average concentrations of U and Th. We ignore 
any contribution from the core. The distribution of radioactive elements are supposed to be 
uniform within each region. Although the concentrations of these radioactive elements are 
considered to be much smaller in the mantle than in the crust, the total geo-antineutrino 
flux from the whole mantle can be comparable to that coming from the crust because of the 
much larger volume. 

Since the total mass of U in the crust is estimated to be Mc(U) = 0.4 x 10^'' kg, and 
Si represents about 15% of the Earth mass, = 5.97 • 10^"^ kg, the other masses in the 
crust and in the mantle can be obtained for each model from the mass ratios they provide. 
In this paper we will fix Mc(U) in the crust to the above estimated value, but we should 
point out that this number is known with 50% uncertainty. The average concentration of U 
in the continental crust is estimated to be 1.7 ppm [12] whereas that in the oceanic one is 
estimated to be 0.1 ppm [13]. We will assume here these numbers as reference values but as 
we will discuss later it is important to know the local variations of U concentration as well 
as as the of the crustal thickness at each experimental site. 

The models can be thus classified by the amount of U they predict for the mantle Mni(U). 
In this work, we follow the classification of models considered in Ref. [6, 7], whose charac- 
teristic will be described below. 

A. Chondritic Earth Model 

The Chondritic Earth model assumes for the Earth's gross composition that of the oldest 
meteorites, the carbonaceous chondrites. The mass ratios for these meteorites [14] are: 
M(Th)/M(U) = 3.8, M(K)/M(U) = 7x10^ and M(U)/M(Si) = 7.3x10"^ [15]. Radiogenic 
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production in the chondritic model easily accounts for 75% of the observed heat flow, about 
30 TW. U and Th provide comparable contributions, each a factor of two below that of K. 
In this model, by taking into account that the total mass of the mantle is 4.1 xlO^^ kg (68 % 
of the total Earth mass), concentration of U in the mantle is about 0.006 ppm. 

B. Bulk Silicate Earth (BSE) Model 

The Bulk Silicate Earth (BSE) model provides a description of geological evidence co- 
herent with geochemical information. It describes the primordial mantle, prior to crust 
separation. The mass ratios here are: M(Th)/M(U) = 3.8, M(K)/M(U) = 10^ and 
M(U)/M(Si) = 9.4 X 10~^. In this BSE model the present radiogenic production, mainly 
from U and Th, accounts for about one half of the total heat flow, 20 TW. The antineutrino 
luminosities from U and Th are rescaled by a factor 1.3 whereas K, although reduced by a 
factor of 5, is still the principal antineutrino source. Concentration of U in the mantle is 
0.01 ppm [16]. 

C. Fully Radiogenic (FR I) Model 

One can conceive a model where heat production is fully radiogenic, with K/U fixed at the 
terrestrial value and Th/U at the chondritic value, which seems consistent with terrestrial 
observations [7] . All the abundances are rescaled so as to provide the full 40 TW heat flow. 
All particle production rates are correspondingly rescaled by a factor of two with respect to 
the predictions of the BSE model. The concentration of U in the mantle for this model is 
about 0.03 ppm. 

D. Modified Fully Radiogenic (FR II) Model 

This model is similar to the previous one, the abundances of U and Th are also rescaled 
with respect to BSE but it assumes that as an extreme case, the total heat flow of 40 TW are 
produced in the Earth only by U and Th, completely ignoring K, as considered in Ref. [6]. 
This is an extreme, geo-antineutrino fluxes much larger than this limit would need serious 
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alteration of source distribution. The concentration of U in tiie mantle for this model is 
about 0.04 ppm. 



III. ANALYSIS METHOD 



Here we explain how we calculate the expected number of geo-antineutrino events, for 
each one of the four models presented in the previous section, at a certain detector position 
for a given exposure. To try to distinguish models we have used a function minimization 
which is explained at the end of this section. 



A. Calculation of the antineutrino fluxes from Th and U decay chains 

We would like to calculate the flux of antineutrino produced in the Earth by the decay of 
a certain isotope that reach a detector for all four geophysical composition models presented 
in the previous section. Following Refs. [7, 8], differential flux of antineutrinos produced in 

— * 

the decay chain of radioactive isotope X that will be measured at a detector position R on 
the Earth, can be expressed by the following integral performed over the Earth volume V^, 

where p(r) is the matter density, C{X,f), tx, mx and nx are, respectively, the concen- 
tration, life time, atomic mass and the number of antineutrinos emitted per decay chain 
corresponding to element X = U, Th. fx{EpJ is the normalized spectral function for ele- 
ment X [17], Pi>^{Et)^, \ f—R\) is the Pg survival probability, which can be averaged out, as 
a good approximation, and bring out from the integral the term: 

(PpJ~l-^sin22eo~0.58, (2) 

where ^0 is the mixing angle responsible for the solar neutrino problem, which is fixed to 
the current best fitted value, sin^ 29q — 0.83, obtained by combining the solar neutrino and 
KamLAND data [18], except in Sec. IV B where ^0 is treated as a free parameter to be 
fitted. We note that the Earth matter effect on geo-antineutrino oscillation is very small 
and therefore, can be safely neglected in our analysis. 

6 



The integration in Eq. (1) can be approximately divided into tliree distinct contributions, 
from continental crust (cc), oceanic crust (oc) and mantle (m), assuming that the matter 
density p(r) and the concentration of X, C{X,r) are approximately constant within these 
three regions and ignoring any contribution from the core, as assumed in Ref. [7], as follows, 

$Pe(^) ^ ^ [m^)icc + mx%, + L-(x)/^] , (3) 

where ~ 6370 km is the radius of the Earth, and we have defined the neutrino luminosity, 
Lp^ , and a factor which depend on the crust thickness distribution over the Earth as well as 

the detector location, /j (i = cc, co, m), as 

li.(X) = ^^i^Wft. (4) 

Txmx 



, / d(/) / dOsinO dr 3, (5) 

Vi Jo Jo JR(s-h{rfi,<i>) P _ d 
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where Vi, Ci{X) and pi are, respectively, the volume, the average concentration of U or Th 
and matter density of region i. We can also write the observable luminosities Ll^{X), in 
units of 10^^ particles per second, for masses in units of 10^^ kg, as 

4^(U) = 7.4M,(U), 4^(Th) = 1.6M,(Th). (6) 

In Eq. (5), h{r,6,(f)) indicates the thickness of the shell (mantle or crust) at the position 
{r,6,(j)). For the mantle we have simply used h{r,6,(f)) = -R©/2, giving = 1.6. For the 
crust calculation (typical thickness being ~ 35 km for continental crust and ~ 7 km for 
oceanic one), we have used the Earth Crust 2° x 2° Thickness map [19], which was obtained 
based on seismology, to compute the /cc,oc at a given detector position. Earth's crust is 
divided into 16200 cells of 2° x 2° where within each cell, the crust thickness is assumed 
to be constant. See Fig. 1 for a schematic illustration of the cell. In Fig. 2 we show some 
iso-contours of Earth crust thickness in the (f) — 6 plane (corresponding to longitude and 
latitude) based on this crustal map model [19] which we will use in this work. 

We observe that in our model assumption the predicted ratio of the geo-antineutrino 
fluxes from U and Th does not depend on details of the integration in Eq. (1) but is given 
by the following simple formula for any model we consider, 

$p,(U) M(U)nuru^ 



$p,(Th) M(Th) nxh TTh' 
7 



1.2 . (7) 



Note that wc arc assuming uniform U/Th distribution so that the ratio of the total amount 
of Th and U, M(Th)/M(U), is considered to be 3.8 in any of the input models under 
investigation. 

We have thus the basic equations for determining radiogenic heat production and neutrino 
flows from models of the Earth composition. In order to have some feehng about the local 
variation of geo-antineutrino fluxes due to the variable crustal thickness, in Fig. 3, we present 
the normalized cumulative geo-antineutrino flux coming from the continental as well as 
oceanic crust (without contributions from the mantle) at Kamioka, Gran Sasso, Homestake 
and Sudbury, as a function of the distance (L) from the source to the detector, computed 
using the information from the 2° x 2° Earth Crust thickness map [19]. We note that due 
to the fact that the abundance of U/Th is much larger in the continental crust than in the 
oceanic one, the geo-antineutrino flux from the former is dominant. 

Since the typical size of the 2° x 2° crustal cell at these detector sites is of the order of 
100 km, we can only reliably compute for distances equal or larger than this distance. In 
Fig. 3 we have extrapolated each curve down to 10 km, as a first approximation. For the 
sake of comparison, we have also plotted the hypothetical case where the entire Earth crust 
has a uniform thickness of 30 km. Neutrino oscillations were not taken into account in the 
calculation of the upper five curves but were included in the lower five. These latter curves 
have been normalized with respect to the no oscillation case. 

From this plot, as far as the geo-antineutrino flux coming from the Earth crust is con- 
cerned, we can see that about 30-40 % of the total crustal antineutrino flux comes from a 
distance within 100 km and about 50-60 % from a distance within 500 km from the detector. 
We note that this flux is about 80%, 65%, 40% and 33% of the total flux, respectively, for 
Chondritic, BSE, FR I and FR II models. This implies that it is very important to know 
rather well, with better than 2° x 2° resolution, the variation of the crustal thickness near 
the detector as well as the local variation of the concentration of radioactive elements. 
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B. Number of geo-antineutrino events 



The number of geo-antineutrino induced events from the decay chain of element X in the 
i-th energy bin, Ni{X), is: 

N,iX) = N,t jdE,^ eiE.J a(^,J , (8) 

where Np is the number of free protons in the fiducial volume of the detector, t is the exposure 
time, e is the detection efficiency, which is assumed to be 100 % for simplicity This integral, 
which is easily computed using the cross section given in Ref. [20], is understood to be 
performed in a certain energy bin. 

We divide the number of events in the positron prompt energy range 0.9-2.6 MeV into 4 
bins with an interval of 0.42 MeV following Ref. [10]. Note that the prompt energy -Eprompt is 
related to neutrino energy as ^'prompt — -Sp^ — (m„— mp)-|-2me = 0.78 MeV where m„, nip 
and rrie are respectively, mass of neutron, proton and electron. The number of events in i— th 
bin is defined as the sum of events coming from U and Th, Ngeo{Ei) = A^i(U) -|- A^j(Th). We 
note that the 1st two lower energy bins contains events coming from both U and Th induced 
geo-antineutrinos whereas the last two higher energy bins contain only events coming from 
Th. 



C. minimization 



We define the function for our geo-antineutrino analysis as follows. 



i=l 



7Vobs(^.) + (0.06 7Vob^(^.))' 

where A'"°^''(£'j), N^^°{Ei) and N^^^°{Ei) are respectively, the number of (simulated) ob- 
servable events (geo + reactor antineutrinos), the number of theoretically expected events 
coming from geo-antineutrinos and from nuclear reactors in the neighborhood of the de- 
tector. As mentioned before, geo-antineutrinos will only contribute to the first four bins. 
We have taken into account in our calculation the total statistical error as well as a 6% 
systematic error only for reactor antineutrino events [10]. Since currently no detector has 
enough data in the energy region interesting for geo-antineutrino observations, we simulate 



N°^^{Ei) according to the Earth composition models, taking into account the reactor an- 
tineutrino background for a given exposure and site. Here, for simplicity, we ignore the 
systematic error for geo-antineutrinos as it is expected to be not so important compared to 
the statistical one for the exposure we consider in this work. 

The above Xgeo '^^^^ be minimized with respect to: (i) $p^(Th) and $p^(U), the total 
geo-antineutrino flux coming from Th and U at a given detector site; (ii) (U) assuming 
Th/U ratio fixed; (iii) sin^2^Q and (U), assuming Th/U ratio fixed. In the first two 
cases the neutrino mixing parameters are fixed to their best fitted values, in the latter only 
AttIq is fixed (see the discussion in the following sections for further details). In addition, 
for some of our analyses, we have also added the function for solar neutrinos, xloh which 
was obtained in Ref. [18]. 

IV. RESULTS 

A. Determination of Geo-antineutrino fluxes 

We first discuss the determination of geo-antineutrino fiuxes at Kamioka as well as Gran 
Sasso sites where the former (latter) has larger (significantly smaller) reactor antineutrino 
background. In this work, we have assumed the reactor antineutrino flux at Gran Sasso 5 
times smaller than at Kamioka site, however, it can actually be even smaller [6]. In this 
subsection, we assume that solar neutrino mixing parameter will be determined with a good 
precision in the future. See, for instance, Ref. [21] for a discussion on the perspectives of 
future determination of the solar neutrino mixing parameters. Due to the effect of oscillation, 
the observable flux suffers a reduction of {Ppe) = 1 ~ |sin^2^Q ~ 0.58, which implies that 
the original flux must be understood as ~ 1.7 times larger than we assumed at the detector 
site. 

Strictly speaking, experiments can only measure the total flux of geo-antineutrinos for a 
certain energy range, regardless of their origin (mantle, crust). They cannot directly access 
the amount of U and/or Th in the crust and mantle separately. Therefore, to be more 
conservative, we first try to consider the total geo-antineutrino fluxes from U and Th as 
free parameters to be fltted by the data for a given input model assumption. Since we still 
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TABLE I: Required exposure in units of kt-yr to identify geo-antineutrinos at the 3 a level, at 
Kamioka and Gran Sasso for each input model we consider in this work. 



Site 


Chondritic BSE FR I FR II 


Kamioka 
Gran Sasso 


1.7 1.3 0.67 0.46 
0.89 0.71 0.38 0.26 



do not have enough data in the geo-antineutrino energy range {Ej^^ < 3.4 MeV), currently 
KamLAND has reported data only for about 0.16 kt-yr exposure [10], we have simulated 
future data according to each one of the Earth composition models. We have used these 
simulated data points as an input in our analysis. By performing a analysis we have 
investigated if the experiments can correctly reproduce the input data and distinguish among 
different models. 

First we have estimated the required exposure, in units of kt-yr to identify the presence 
of geo-antineutrino flux for an assumed model input, leaving the two components of the 
geo-antineutrino flux, coming from U and Th, completely free in our fit. The results are 
shown in Table I. The numbers in the table are the detector exposure, at each site, required 
to rule out or accept each model at 3 o" level. The actual fiducial volume of KamLAND 
(Borexino) is about 0.4 (0.1) kt, which means one has to multiply by a factor 2.5 (10) the 
number found in Table I in order to translate it to the actual detector exposure time. As 
naturally expected, the larger the input flux, the easier to identify the model, and therefore 
we need a smaller exposure to rule it out. If, for instance, a detector such as KamLAND, 
does not see geo-antineutrinos at the 3 a level after 0.46 kt-yr exposure, it can rule out 
the FR II model, but at this point it cannot say anything about the other models. On the 
other hand, if KamLAND measures geo-antineutrinos at this exposure, then FR II has to 
be interpreted as the preferred model, and so forth. 

In Fig. 4, we present in the $Pg(U) — $p^(Th) plane, to which extent we can determine 
geo-antineutrino fluxes in the presence of neutrino oscillations for given input model fluxes 
and exposure (in units of kt-yr) at the Kamioka site. In all cases we are able to correctly 
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reproduce the input fluxes as the best fit points (indicated by stars) in our fit, which, 
however, does not necessarily mean we can distinguish models with enough significance (see 
the text below). Unfortunately it seems to be extremely difficult to distinguish between 
BSE and Chondritic models, independently of the detector exposure we have considered 
here. Prom this plot, we can conclude that for 1 kt-yr of exposure, it is not possible to 
distinguish among the four models considered. However, after 3 kt-yr of exposure, one start 
to have some sensitivity to distinguish models, namely, Chondritic/BSE from FR II. After 
the maximal exposure considered, 6 kt-yr of data, one can almost always distinguish FR II 
from any of the other models. It is worthwhile to note that even for the maximal exposure 
and maximal flux (FR II) considered here, it is not possible to exclude a null Th induced 
geo-antineutrino flux, whereas it is sometimes possible to exclude a null U induced one. 

In Fig. 5, we present the same plot but for the Gran Sasso site, assuming a significantly less 
reactor neutrino background, 5 times smaller than at the Kamioka site. In this case, one can 
achieve the same sensitivity for the flux determination with much less exposure compared to 
the Kamioka site. Roughly speaking, x kt-yr exposure at Gran Sasso correspond to (2 — 3)a; 
kt-yr at Kamioka, which we can see by comparing Figs. 4 and 5. The general behavior of 
the allowed regions are very similar to that for Kamioka site apart from this difference. 

In Fig. 6, we show with which precision the total U induced antineutrino flux can be 
determined by the experiments, further imposing that the ratio between the total U and Th 
mass in the Earth is fixed to be constant, M(Th)/M(U) = 3.8, in the fit. Note that this 
reasonable assumption, based on meteorite and Earth surface data, is not very important in 
distinguishing models but allows for a better precision in flux determination. In this case, 
we have only one free parameter, $p^(U), to be fitted. In 6 kt-yr at Kamioka, the total U 
flux can be determined within less than 10 % uncertainty, independently of the model. The 
precision it can be determined at Gran Sasso, after 2 kt-yr exposure, is a bit over 10 %. 

Let us now translate the fiux uncertainties we have estimated into heat uncertainties, 
in order to compare our results with those of Ref. [9]. As one can easily understand from 
Eq. (8) of that paper, we can write the total heat, H, as a. function of the geo-antineutrino 
flux as 

H = a$p, - b , (10) 
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where a and b are positive constants satisfying a$p^ > b. This implies that the relative heat 
uncertainty is always greater than the relative flux uncertainty, that is, 



AH 




(11) 



H 



> 



We have roughly estimated that the flux uncertainties of the order of 10 %, that can 
be achieved at Kamioka with 6 kt-yr exposure, correspond to 15-30% in heat uncertainties, 
consistent with the results presented in Ref. [9]. 

B. Dependence on the solar mixing angle ^0 

So far, we have flxed the solar mixing angle ^0 to its current best fltted value. How 
our ignorance about the precise value of the solar mixing angle can aggravate our results? 
To answer this question, we perform a flt leaving sin^ 26q as a free parameter. For this 
purpose we will combine present solar neutrino data with future simulated reactor and geo- 
antineutrino data. 

As a flrst step, to see how the observation of geo-antineutrinos can affect the determina- 
tion of the mixing angle, we have analyzed the present KamLAND reactor antineutrino data 
(17 bins) allowing for geo-antineutrino contributions in the flrst 4 bins, as performed by the 
authors of Ref. [8] leaving Th and U antineutrino fluxes free but imposing that the ratio 
between the total U and Th mass in the Earth is fixed to be constant, M(Th)/M(U) = 3.8, 
as in the end of the previous section. In Fig. 7, we show the result we have obtained. In 
agreement with the result in Ref. [8] , the allowed region becomes somewhat smaller when we 
include events which can be interpreted as geo-antineutrinos. However, we have confirmed 
that if U and Th contributions were treated as independent free parameters, the allowed 
region would not shrink as much, as recently pointed out by Inoue in Ref. [10]. Some events 
observed in the energy range -Eprompt < 2.6 MeV can be attributed to geo-antineutrinos, but 
the claim that geo-antineutrinos have been observed can not be made at this point due to 
small statistics (see discussion in previous section). 

We further proceed, by combining the KamLAND data with the current solar neutrino 
one. The final result is shown in the right panel of Fig. 8, where we have also presented the 
result without the geo-antineutrino constraint in the left panel taken from Ref. [18]. The 
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allowed region becomes somewhat smaller but essentially has not changed. This is quite 
understandable as solar neutrino data are dominating in the determination of the mixing 
angle at this point. 

Finally, we combine the simulated future reactor and geo-antineutrino data with the 
present solar neutrino data. In order to see how the discrimination power of composition 
models depends on our ignorance of the precise value of the solar mixing angle, for a given 
input value (fixed to be the current best fit value) as well as $p^(U), we perform a fit 
leaving $p^(U) and ^0 as free parameters. In this fit, we assume the ratio between the 
total U and Th mass in the Earth to be constant, M(Th)/M(U) — 3.8, for simplicity. This 
is quite reasonable for our purposes here as the discrimination power among models does 
not seem to depend very much on this assumption, as we have discussed in the previous 
subsection. 

In Figs. 9 and 10 we plot the allowed region for "I>Pg(U) and sin^ 2^q for the Kamioka and 
Gran Sasso sites. The reduction of the allowed range of sin^ 2Qq from 1 to 6 kt-yr exposure 
is essentially due to reactor antineutrinos. We observe that our ignorance on the exact value 
of ^0 does not influence very much the determination of $pg(U), so our conclusions in the 
previous subsection do not essentially change. 

V. DISCUSSION AND CONCLUSION 

The amount of radioactive elements in the Earth is not well known. The total quantity of 
U and Th in the Earth, however, can be directly measured by neutrino detectors. Presently 
KamLAND data imply that the Earth radiogenic heat output can be anything between 
and 110 TW [10] and all Earth composition models are compatible [8]. We have studied 
how this can be improved by future data. 

We have investigated to which extent U and Th geo-antineutrino fluxes can be determined 
by neutrino detectors in a decade of exposure. We have considered KamLAND at Kamioka 
and Borexino at Gran Sasso as our reference detectors and sites. We have found, as we 
showed in Table I, that within a few years with a relatively small amount of exposure, 
it is possible to establish the presence of geo- antineutrinos unless their flux is signiflcantly 
smaller than expected. However, to discriminate among different Earth composition models. 
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considerably longer exposure is required. We found that in 6 kt-yr at Kamioka, the total 
U flux can be determined within less than 10 % uncertainty, independently of the model 
we considered. The precision it can be determined at Gran Sasso, after 2 kt-yr exposure, 
is shghtly larger than 10 %. We note that at Gran Sasso the same sensitivity to geo- 
antineutrino flux determination can be achieved with substantially smaller exposure due to 
much lower reactor antineutrino background. 

We observe that our ignorance on the exact value of ^0 does not influence very much the 
determination of the geo-antineutrino flux. However, it is very important to know the local 
variation of the Earth crustal thickness as well as concentration of U and Th in the region 
close to the detector with better than 2° x 2° resolution, to be able to accurately translate 
the measured flux into amounts of U and Th. 

The determination of the radiogenic component of the Earth heat generation is of great 
geophysical interest. Experiments such as KamLAND and Borexino can open a new window 
to survey the internal structure and dynamics of our planet, leading to the birth of neutrino 
geophysics. 
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Noth Pole 




FIG. 1: Schematic illustration of the k-th cell of the Earth crust in which its thickness is assumed 
to be constant. 



Earth Crust Thickness Map 




FIG. 2: Iso-contours of the Earth crust thickness based on the global 2° x 2° crustal map [19] we 
adopted in this work. 
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10 10^ 10^ 10^^ 

L (km) 

FIG. 3: Normalized cumulative geo-antineutrino flux, with (lower five curves) and without (upper 
five curves) neutrino oscillation, coming from the continental as well as oceanic crust (without 
contributions from the Mantle) at various positions on the surface of the Earth as a function of the 
distance (L) from the source to the detector. These curves were computed using the information 
from the 2° x 2° Earth Crust thickness map [19] . For the sake of comparison, we have also plotted 
the hypothetical case where the entire Earth crust has a uniform thickness of 30 km. 
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FIG. 4: Region allowed at 1, 2 and 3 a in the $Pg(U) — $Pg(Th) plane, for 4 different models (in 

first 4 columns) for 1, 3 and 6 kt- yr of exposures (in different rows). In the 5th column, we have 

superimposed 3 a allowed regions of 4 different models. In the second row, for each model, the 

input values for the geo-antineutrino fluxes computed according to Eq. (1) using the information 

from the Earth Crust 2° x 2° Thickness map [19] are indicated. Both of <I>p^(U) and <I>p^(Th) fluxes 

are allowed to vary freely. Dashed lines indicate the expected U and Th induced antineurino fluxes 

if the Th/U mass ratio is assumed to be constant, M(Th)/M(U) = 3.8. The input points, which 

practically coincide with the best fit, are indicated by stars. The input fluxes are estimated in the 

presence of neutrino oscillations with Am^ = 7 x 10~^ eV^ and sin^ 2Q = 0.83. 
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FIG. 5: Same as in Fig. 4 but for the Gran Sasso site. 
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FIG. 6: The allowed ranges of geo-antineutrino flux from U are plotted for given different input 
flux assuming M(Th)/M(U) = 3.8 for KamLAND for 1 and 6 kt- yr of exposures and for Borexino 
for 0.5 and 2 kt- yr exposures. 
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FIG. 7: Regions in (sin^ 2^q, Am^) plane allowed by KamLAND data alone (without solar neutrino 
data) for different thresholds ii^prompt > 0.9 MeV (including geo-antineutrino candidates) and 
-^prompt > 2.6 MeV (without geo-antineutrinos contribution). We have imposed M(Th)/M(U) = 



3.8 in the fit. .|q ^ 



KL without + Solar Data 



KL with z/,e„ + Solar Data 




FIG. 8: Region allowed by all the solar neutrino experiments combined with KamLAND data with 
different thresholds E'prompt > 2.6 MeV (left panel) and ii^prompt > 0.9 MeV (right panel). We have 
imposed M(Th)/M(U) = 3.8 in the fit. 
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KamLAND + Solar Data combined 
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FIG. 9: Region allowed at 1, 2 and 3 a in the ^p^{\J) — sin^ 29q plane, for 4 different models (in 
first 4 columns) for 1, 3 and 6 kt- yr of exposures (in different rows). In the 5th column, we have 
superimposed 3 a allowed regions of 4 different models. Here the flux ratio of (U) and <I>p^ (Th) 
is fixed to be constant. 
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Borexino + Solar Data combined 
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FIG. 10: Same as in Fig. 9 but for the Gran Sasso site. 
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